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Abstract. Lo, Sargent, and Young (1993) have recently published an anal- 
ysis of the HI kinematics of nine faint dwarf galaxies. Among other things, 
they conclude that the masses of these systems, as deduced by the modified 
dynamics (MOND) from the observed velocity dispersions, are systematically 
smaller than even the HI masses that are observed in these systems, by a 
factor of ten or more. Such a state of things would speak strongly against 
MOND. We show here that the MOND mass estimator used by Lo et al. is 
smaller than the proper expression, by a factor of about twenty. We derive 
the proper mass estimator as an exact virial-like relation between the 3-D rms 
velocity, (i?^), and the total mass, M, of an arbitrary, self-gravitating system, 
made of light constituents, that is everywhere in the very-low-acceleration 
regime of MOND. This reads M = {9/A){v'^y /Gao- (For a system that is 
not stationary, (f^) involves also an average over time.) We do this in the 
Bekenstein-Milgrom formulation of MOND as a modification of gravity. This 
relation has been known before for the special case of a stationary, spherical 
system. We further generalize this relation to cases with constituent masses 
that are not small compared with that of the whole system. We discuss vari- 
ous applications of the M — v relation; inter alia, we derive an expression for 
the two-body force law in the large-distance limit. With the correct estimator 
the predictions of MOND are, by and large, in good agreement with the total 
observed masses (the observed gas mass plus a stellar mass corresponding to 
an M/L of order one solar unit). 



1. Introduction 

Lo, Sargent, and Young (1993) (hereafter LSY) have recently published 
the results of HI observations of nine intrinsically faint dwarf galaxies. Among 
other things, they determine total HI masses. They also estimate and discuss 
the total masses implied by the modified dynamics (MOND). They find that 
the latter are systematically smaller than the HI masses "by a factor of ten or 
more" -an unacceptable state of affairs (for the theory). 

LSY have used a MOND mass estimator that is a mistaken adaptation 



of the relation 



M = V^/Gao, 



(1) 



between the total mass, M, of a body, and the asymptotic rotational velocity of 
test particles around it, (Milgrom 1983a,b); Gq is the acceleration constant 
of the theory. 

Equation (1) is exact in MOND, but LSY have used an analogous relation 
with Voo replaced by Vh, the half width at half maximum of the integrated HI 
line profile-which is a measure of the one- dimensional rms velocity for the 
whole galaxy. LSY thus use 

M, = vt/Gao, (2) 

We show below that gives a far underestimate of the mass. We derive 
a general mass- velocity relation that holds in the non-relativistic, Bekenstein- 
Milgrom (1984) formulation of MOND as a modification of gravity; this rela- 
tion is 

M=^{v'y/Gao. (3) 

Here, (v^) is the three-dimensional mean-square velocity of the whole system 
(averaged over time if the system is not stationary), and M is its total mass. 
Relation (3) applies for an arbitrary, self-gravitating system, whereby the 
acceleration is much smaller than Oq everywhere. Beside the factor of 9/4 
missing in expression (2) , a large factor of disparity (9 in the isotropic case) 
is present because of the difference between 3-D and 1-D velocities, leading 
to an undrestimate of the masss by a factor of about twenty, in the isotropic 
case. 

Gerhard and Spergel (1992) have derived relation (3) for the special 
case of a stationary, spherical system; it was derived, for the yet more special 
case of a sphere with constant radial and tangential velocity dispersions, by 
Milgrom (1984). 



II. The MOND mass estimator 

Consider a self-gravitating system that is composed of various particle 
species with masses and distribution functions fk{f,v,t). As usual, take 
the time derivative of the quantity 

5: /.^r (4) 
k 



(itself, the time derivative of the trace of the moment-of- inertia tensor). 




(By Liouville's theorem the time derivative of a quantity of the form 
J d^rd^vfq{f,v,t) is J d^rd^vfq.) The first term in eq. (5) is the momen- 
tary, mass- weighted, 3-D, mean-square velocity, {v^), multiplied by the total 
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mass, M. In the second term we put a = —Vcf) (here it is assumed that gravity 
is the only important force), where (f) is the (MOND) gravitational potential, 
to obtain 

Q = M{v^){t) [d'r d'v rrik h{r,v,t) f • V0. (6) 

k 



As a — — V0 depends only on r, the v integration, and the sum over species 
can now be performed to give the standard result 

Q = M{v^){t) - J d^r p{f,t)f- V0. (7) 

We now specialize to MOND, and here enters the assumption of self gravity: 
the density p is the source of the gravitational potential. We substitute p(r, t) 
from the field equation of Bekenstein and Milgrom (1984), 

V • [p{x)V(f)] = AnGp, X = |V0l/ao, (8) 
in eq. (7) ; then, integrating by parts, we get 

Q = M{v^){t)-^ j^li{x)T-Vci>Vcj>-ds+^ j pi{x)Vcj>-V{r-V<j>)d'T. (9) 

The first integral is over any surface, S, encompassing all the mass; we take 

it at infinity. As r goes to infinity V(/> becomes {MGaoY^'^r/r'^ (see Beken- 
stein and Milgrom 1984), and p{x) becomes x, so the second term becomes 
—M{MGaoY^^- The second integral can be calculated thus 

/2 = j ii{x)V^-V{r-V^) d^r = j pi{x)V ^ ■ [{r ■ V)V ^ + {V ^ ■ V)r\ d^r = 

J ^pix)f- V[{V<pf] + ;u(x)(V0)' d^r. (10) 

Let F{y) be such that p{x) = F'{y)\y^^2 [F(x^) is the Lagrangian density 

for the potential (see Bekenstein and Milgrom 1984)]. The first term in the 
integrand of eq. (10) can then be written as 

„2 2 o 2 

|f.VF=|v.(Fr)-^F. (11) 

The first divergence term may again be written as a surface integral at infinity 
that contributes (l/3)M(MGao)^/^ to the right-hand side of eq. (9) . We thus 
end up with the relation 

Q/M = {v')it) - liMGaoY^' - J ^F{x') - p{x)x' d\ (12) 



(with X = \V(f)\/ao). Take now the long-time average of eq. (12) . That of the 
left-hand side vanishes, as Q is finite at both ends of time, and we get 

F) = ^(MGao)^/^ + ^^ lpW)-l&^d\ (13) 

where an overline signifies the long-time average. This relation is exact and 
general (i.e independent on how close to the MOND regime we are). In the 
Newtonian limit (tto goes to 0) it gives the usual virial theorem. In the very 
limit where the accelerations are always and everywhere in the system much 
smaller than Oo, the integrands in eqs. (12) (13) vanish [because in the limit 
a; — > we have F{x^) {2/3)x^, and iJ,{x) x], and we obtain 

W) = l{MGaoy/'. (14) 

When the system is stationary (and sometimes in more general cases-such as 
systems that are stationary in some rotating frame) the momentary value of 
(f^) remains constant and it can then be used in eq. (14) . As we cannot 
determine time averages for relevant astrophysical systems, we usually make 
the simplifying assumption that the system is stationary and use relation (14) 
with the observed momentary value of (f^). 

Equation (14) is the trace of the following relation , which holds under 
the same conditions, and which can be derived along the same lines: 

= liMGaoy/%j + J '{S~^3^~^~JW^ d'r. (15) 

There are possible formulations of MOND other than that of Bekenstein 
and Milgrom (1984). The present M — v relation is probably not exact in all 
of them. It is, however, a good indicative estimator, and, in any case, the best 
we have at the moment. We note, for example, that in the class of MOND 
theories based on modification of the law of inertia, discussed by Milgrom 
(1993), the M — V relation is still exact for stationary, spherical systems, 
whose constituents move on circular trajectories (in the deep MOND limit). 

In many instances, what we would like to use as test particles are them- 
selves sub-systems with internal structure, or with masses that are not small 
compared with M. We must then be careful in treating the subsystems. Let m 
be the mass of a body we want to treat as a structureless constituent, and let 
I be its size. The assumptions underlying our derivation of the M — v relation 
seem to require the following: (i) The constituent itself is kept together by 
gravitational forces balancing internal motions (say of some perfect fluid), (ii) 

In order to have |V(/)| <S Oo everywhere (within the body as well as outside) 
we must have mG/P ^ ao- (iii) If L is the length scale over which the field 
varies appreciably, in the vicinity of the body, we must have £ -C L. (iv) As we 
want to use only the centre-of-mass velocity of the body in calculating (v^), 
and neglect the internal velocities, we must have m ^ M . 

This set of requirements is too restrictive for the M — v relation to be 
of much application: Condition ii, for instance, would bar stars from being 
legitimate constituents; condition i would bar atoms and elementary particles, 
etc. (condition iii is quite benign). Fortunately, it is possible to weaken 
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conditions i,ii greatly, or, indeed, to eliminate them altogether. Also, the 
M — V relation may be generalized to circumvent condition iv. 

We note first that when all the would-be constituents have m -C M 
requirements i,ii may be disregarded. To see this, replace all bodies that do 
not satisfy i and/or ii by ones that do; i.e., by ones having the same mass, but 
a size satisfying mG/£l <^ Gq, and that are made of a perfect fluid held by 
gravity. Because m -C M, can be chosen small enough so that <^ L. The 
new system does satisfy all the assumptions, and the M — v relation holds for 
it. However, If m is small enough there exist a length a such that a -C L on one 
hand, and such that at a distance a from m the contribution of the latter to 
V(/) is small. Under these conditions, Bekenstein and Milgrom (1984) showed 
that the centre-of-mass motion of the body is oblivient to its internal structure 
and mass, and it behaves like a test particle. The centre-of-mass motions of 
our replacement bodies are, then, like that of the original. The correction on 
this very-small-mass case is of order (m/M)^/^. In addition, we want to be 
able to neglect the internal velocities in the replacement constituent bodies; 
again this is permitted because m <C M, and the higher order correction is 
here also of order (m/M)^/^. It is thus valid to apply the M — v relation to 
the original system without having to worry about the internal makeup of the 
constituents. 

When at least some of the sub-system's masses are not very small com- 
pared with M, the M — v relation, as given by eq. (14) is not valid. Given 
the masses and sizes of the constituents, we consider the limit where the dis- 
tances between them are very large-in keeping with our working in the extreme 
MOND limit. We described before, replace the masses with "live" , self- 

gravitating ones made of a perfect fluid. Relation (14) can now be applied 
but we must reckon with the internal velocities of the "live" replacements, 
which does not necessarily have a counterpart in the original constituent. In 
the limit of large distances, each replacement body itself satisfles eq. (14) , 
as can be shown. (We assume that the internal acceleration produced by the 
constituent dominates over that of the rest, at a distance where the asymp- 
totic form of the isolated-mass potential-as described below eq. (9) -already 
holds. This is valid in the limit that we consider here.) The prescription is 
thus as follows: if is the mass of the ith particle, and Vi its centre-of-mass 
velocity, then the mean square velocity in the M — v relation is to be taken as 



where the second term in parentheses is the internal rms velocity within mj. 
This leads to the relation 



where (v^) now includes only centre-of-mass velocities. 

To demonstrate the neccesity of this crucial modiflcation, consider a sys- 
tem made of two masses m and M ^ m, in a circular orbit around each other. 
A blind application of relation (14) will give MGao ~ (9/4) (m/M)^^'*, where 
V is the speed of the smaller mass. This, however, is quite wrong: MOND tells 
us that, in fact, MGao ~ v"^- The reason for the failure of eq (14) , in this 




(16) 




(17) 
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case, is that most of the contribution to (f ^) should have come from intrinsic 
motions within M, which are neglected when we treat it as a test particle. If 
we apply eq. (17) instead, we get the correct result (as the first order term in 
m/M). 

The relative correction introduced by eq. (17) is bounded by (m^/M)^/^, 
where rrix is the maximum mass of a constituent. When all the masses are 
equal the relative correction is (m/M)^/^. 

At any rate, when the accelerations internal to the subsystems are small 
compared with ao, it is best not to consider them as test particles, but include 
in (f ^) the full rms velocity for the system. 

III. Some applications of the mass-velocity relation 

First, and perhaps foremost, the general M — v relation puts on firmer 
and wider footings the MOND prediction of a Faber- Jackson like relation for 
all self-gravitating low-acceleration objects. This includes disk galaxies, as we 
saw; for these, the M — v relation comes in addition to the (differnt) relation 
between total mass and asymptotic rotational velocity (1) . These relations 
tell us that we can use in the TuUy-Fisher relation either the rms roatational 
velocity (for low-a galaxies) or the asymptotic, rotational velocity (for all 
disks) . We then expect relations with the same slope, but different intercepts. 
In fact, almost all known astronomical objects, on the scale of galaxies and up, 
have internal, mean accelerations that are about ao or smaller, and the M — v 
relation should approximately hold universally. In applying it to observed 
objects we note again the uncertainties in the deduction of the 3-D velocities, 
and also the error introduced if we assume that a system is already stationary 
when, in fact, they are not. (and take its momentary value of (v^) to be the 
time-average value). 

The M — V relations affords significant shortcuts to the derivation of 
the MOND forces acting on bodies in a few configurations. If the gravita- 
tional forces can be balanced by rotating a configuration, rigidly, with some 
frequencey a;, about some known centre, then, uj may be determined from re- 
lation (17) , and from it the forces on all the masses may be gotten. Perhaps 
the most interesting application is a derivation of an analytic expression for 
the two-body force in the Bekenstein-Milgrom formulation, in the long dis- 
tance limit: Consider two masses mi, and m2, at a distance I from each other: 
a system that can be supported by rotation about the centre of mass. In the 
limit of large £ an expression for the force, F, can be conveniently put in the 
form 

F(mi,m„£)= ^ (-^^\'\{m,/m,\ (18) 

with 

A{q) = \q-\l + + qf' - q''' - 1]- (19) 

Numerical results for the dimensionless function A{q), for a few values of the 
mass ratio, q, are given by Milgrom (1986), and agree with those determined 
from the analytic expression (19) . We note that A[q) varies rather little across 
its full range, from A=lsXq = {)ioA^ 0.8, q = 1. {A is, of course, 
symmetric under the exchange of the two masses q <-> q~^, because the forces 
on the two masses are equal.) 
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Similarly, we can calculate the force (of a more academic interest) on any 
of n equal masses M/n, on the vertices of a regular polygon, of diameter 2r, 
with a mass m at its centre. 

F{M, m, n,r) = —{M + m)^^^{GaoY^^ 
3n 

This was checked numerically, by R. Brada (private communication) , for some 
configurations. 

It is useful to specialize the M — v relation to rotationally supported disks 
(e.g. a very low acceleration disk galaxies). If S(r) is the surface density, and 
v{r) the rotation curve, then eq. (14) implies the relation 

o POO 

^M^/^{Gaoy/^ = / 27rrS(r)^;'(r) dr, (21) 
3 Jo 

where M = 27rrE(r) dr, is the disk's mass. If we add a point mass m at 
the centre, representing a bulge, say (within a radius that is small compared 
with the dimensions of the disk), eq. (17) tells us that the relation is now 

-[(M + m)^/'-m2/'](Gao)^/' = / 27rrS(r)^;'(r) dr. (22) 
"J Jo 

If we consider a special case of a thin ring of mass M and radius r, with a 
point mass m in the centre, we can use eq. (22) to calculate the gravitational 
force on the ring (M times the force on a unit mass): 

F=l{M + mY/^{Gaoy^^ 

O 

(This is also a special case of eq. (20) , taking there the limit of nF for 
n — > oo.) 



1 - 



^3/2^ ^-1/2^3/2 

(M + m)3/2 



r-\ (20) 



m 



M + m 



3/2 



(23) 



IV. A revistation of the masses of the dwarfs 

The 3-D rms velocity, which enters the M — v relation is not, in itself, 
directly observable, and we have to express (v^) in terms of the observed, 
integrated, 1-D rms velocity, fy. This, as usual, involves assumptions on the 
velocity structure of the system. If is independent of the direction of the 
line of sight-as when the velocity distribution is isotropic-we have (f ^) = 3f y , 
and the mass estimator for this case is 

M = ^(^;/G'a«). (24) 

If the motions in the system are confined to a plane that is at an inclina- 
tion i to the line of sight, and if, in that plane, the overall velocity distribution 
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is isotropic-as is the case for any axisymmetric systems, such as a rotation- 
supported disk-then (v^) = and 

M=^^.{vyGao). (25) 

We see that the LSY estimator could, in the case of flat systems, be even more 
discrepant with the correct one. For example for i < 45° we get a factor of 
discrepancy > 36 (for i < the factor is > 144). 

Two factors conjoin to give the large disparity between the naive mass 
estimator used by LSY, and the correct expression. The estimator given in 
eq (1) , which LSY mimic, uses the asymptotic circular velocity around the 
system. However, the (3-D) rms velocity in a system is smaller then the asymp- 
totic rotational velocity corresponding to system's mass. Also, the measured 
line-of-sight rms velocity is of course smaller then the 3-D one. Since these 
ratios enter in the fourth power, the resulting disparity is large. As v\\ is 
approximately Vh we see that the mass estimator used by LSY differs by a 
large factor from the correct one. [In fact, f y may be somewhat smaller than 
Vft-according to LSY Vh ~ 1.18v||-but the correction introduced by this differ- 
ence is compensated by the fact that LSY use the outdated pristine estimate 
of Go given by Milgrom 1983b, when the best value today, which is based on 
detailed studies of rotation curves, is about a factor of two smaller (Begeman, 
Broeils, and Sanders 1991)]. 

In default of more detailed information on the galaxies studied by LSY, 
the best we can do to estimate their MOND masses is to use relation (24) . 
This means that we simply have to multiply LSY's MOND masses by a factor 
of twenty. Because deduced MOND masses are subject to large errors, as they 
depend on the fourth power of the velocities, we feel it is more appropriate to 
use the "luminous" masses in relation (24) in order to calculate from them the 
velocity dispersions that MOND implies, and compare these with the observed 
line-of-sight velocities v\\ (designated ctt by LSY). We do this in Table 1 where 
we give f ^ = 0"^ along with the total line-of-sight velocity dispersions predicted 
by MOND on the basis of the stationary M—V\\ relation (24) , for the isotropic 

case. These are calculated for three values of the stellar M/L: a = 0.5, 1., 4 
solar units, and are designated v\\{a).{'Wc multiply the HI mass by a factor 

of 1.3 to account for He.) We use the value ap = 1.2 lO^^cm/ (Begeman, 
Broeils, and Sanders 1991). We see that the observed fy falls within the 
predicted range of velocities corresponding to the range of M/L values we use, 
in all but two galaxies, where MOND predict higher velocities than observed 
(the quoted errors in are about ±20%). 

We must note that the comparison for individual cases is subject to large 
uncertainties due to the following factors: (i) We do not know that the systems 
under study are in a stationary state, so that their momentary, observed rms 
velocity equals the long-time average. In fact, LSY state that the velocity fields 
of some of the galaxies suggest radial motions such as expansion or contraction. 
Non-stationary systems spend more time in a state of lower-than-average rms 
velocity. If departure from stationarity is important we would expect that the 
velocities deduced from the M — v relation-representing long-time averages- 
would be, by and large, larger than the momentary, observed ones, (ii) Our 
M — V realtion assumes that the HI is the dominant mass component (i.e. 
is self gravitating). This is clearly not so in some of the cases (for which 
LSY find Mh/Lb < 1); notably in LGS-3, and DD0216, but also in DD069, 
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and DD0155). In these cases, the deduced velocities depend strongly on the 
distributions of both the stellar, and the HI masses which we do not know. 
Here there is an additional uncertainty in the value of the observed masses 
due to that in the stellar M/L, which we try to cover in Table 1. (iii) The 
distances to the galaxies in the sample arc poorly known according to LSY. 
The value of the deduced velocities is proportional to the square root of the 
observed distance, (iv) The observed line-of-sight velocities need not be 3"^^^ 
of the 3-D rms velocities as we assume in applying the M — v\\ relation for the 
isotropic case. 
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Name 




Lb 


^11 


t;|l(.5) 


^11 (1) 


^11 (4) 


LGS-3 


0.02 


0.07 


7.5 ±1.3 


4.7 


5.3 


7.0 


UGC4483 


3.7 


2.3 


13.0 


14.5 


15.4 


18.2 


DD069 


3.6 


3.6 


8.0 


15.0 


16.0 


19.7 


CVn dwA 


8.1 


1.1 


9.5 


17.2 


17.4 


18.5 


DD0155 


0.2 


0.23 


10.5 


7.4 


7.9 


9.9 


DD0187 


5.0 


2.3 


13.0 


15.6 


16.2 


18.7 


Sag DIG 


0.8 


0.25 


8.0 


7.8 


10.1 


11.3 


DDO210 


0.3 


0.13 


8.0 


7.7 


8.0 


9.2 


DD0216 


1.3 


6.9 


10.0 


14.2 


16.1 


21.9 



Table 1. Observed masses, luminosities, and line-of-sight velocity dis- 
persions for the dwarfs, along with velocity dispersions predicted by MOND 
v\i{a), for assumed M/L = o: in solar units. All masses in units of 10''' M©, 
luminosities in units of 10^ Lq, and velocities in km/s. 
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